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We apply the worldline formalism and its numerical Monte-Carlo approach to compu- 
tations of fluctuation induced energy-momentum tensors. For the case of a fluctuat- 
ing Dirichlet scalar, we derive explicit worldline expressions for the components of the 
canonical energy-momentum tensor that are straightforwardly accessible to partly ana- 
lytical and generally numerical evaluation. We present several simple proof-of-principle 
examples, demonstrating that efficient numerical evaluation is possible at low cost. Our 
methods can be applied to an investigation of positive-energy conditions. 
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1. Introduction 

It has been over five decades since Casimir predicted a negative interaction energy 
density and a corresponding attractive force'' between two perfectly conducting 
plates in a vacuum which impose boundary conditions on the quantized electromag- 
netic field. Although the theory of the Casimir effect has become a well developed 
field in recent year^, it still serves as a paradigmatic test case for exploring the 
diversity of phenomena within quantum field theories. In the first place, interac- 
tion energies, pressures and forces caused by boundary conditions (BCs) imposed 
on a quantum field are of particular relevance, as they constitute phenomenological 
observables of the Casimir effect and can be directly related to experiment. From 
a more theoretical perspective, this information is more generally encoded in the 
expectation value of the energy-momentum tensor (EMT) of the fluctuating field. 

Beyond the direct phenomenological application, the energy-momentum tensor 
is also at the center of interest since it serves as a source term for the gravita- 
tional field if the quantum field is coupled to gravity. For instance, the energy- 
momentum tensor is needed in order to study how the Casimir energy falls in a 
gravitational fielcP. For usual Casimir force calculations, the fact that Casimir en- 
ergies can generically be negative is not a problem at all, as it is similar to negative 
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binding energies in bound state systems. However, in connection with general rel- 
ativity, negative energy densities can be puzzling. Such energy densities may be 
in conflict with restrictions imposed on the energy-momentum tensor in order to 
avoid exotic phenomena such as superluminal travel, wormholes or time machines. 
Indeed, the Casimir effect violates particular restrictions such as the weak energy 
condition (WEC), (T^uVf^V^) > with being a timehke vector, or the null en- 
ergy condition (NEC), {f^uV'V) > with V being a nuU vectoJUl By contrast, 
both conditions are obeyed by classical physics. 

A somewhat weaker condition that is still sufhcient to rule out exotic phenom- 
ena is the averaged null energy condition (ANEC) which requires that the null 
energy condition has to hold only when integrated along a complete geodesic. The 
ANEC has been found to be satisfied in all Casimir examples studied so far on 
flat Minkowski space (even though it can be violated on compact flat or curved 
space^. An investigation of the ANEC in specific Casimir configurations is typi- 
cally involved, as certain components of the energy-momentum tensor need to be 
known along a complete geodesic. However, the physically interesting geodesies 
which appear to collect negative energy densities along the geodesic are typically 
directed towards a surface and thus hit the surface at some point. In such a case, 
the average acquires a generically large positive contribution from the surface itself 
which - though nonuniversal - is likely to exceed the negative Casimir contribu- 
tions by far. Relevant configurations that avoid the discussion of large nonuniversal 
contributions therefore need to have holes through which a geodesic can pass. In 
fact, the case of a single plate with a hole has been found to obey the ANEC!^. 

As is obvious from these considerations, a general study of energy conditions in 
Casimir configurations requires a theoretical framework that is capable of dealing 
with the energy-momentum tensor in arbitrary geometries. Standard computations 
of Casimir energy momentum tensor^lHMIIIIl! are usually based on mode sum- 
mation/expansion, image charge methods or similar techniques which are suited 
for specific simple geometries. In the present work, we report on progress using 
the worldline approactfl^ to the Casimir effect for the case of a Dirichlet scalar 
fielcCHill which has meanwhile been used to study a variety of nontrivial Casimir 
configuration^lEIElEH We show that the worldline approach can straightforwardly 
be generalized to study Casimir energy-momentum tensors. We discuss the similar- 
ities and differences with conventional interaction energy computations and verify 
the resulting algorithm with the help of known simple examples. 



2. Worldline Formalism for composite operators 

We investigate the energy-momentum tensor of a minimally coupled quantum scalar 
field with zero mass which shall be defined on a d + 1-dimensional domain V. 
On the boundary dT> that we consider as static we impose Dirichlet boundary 
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conditions. The canonical energy-momentum tensor operator of $ reads 



T^y{x,t) := lim^ 



(1) 



where we have employed a point-sphtting procedure in the spatial coordinates for 
regularization. The background potential a{x) can in principle be arbitrary but 
will later specifically be used to impose Dirichlet BCs. For our study of the energy 
conditions, we need to compute the vacuum expectation value (t^i,'^ of the EMT. 
For instance, for the null energy condition, we have to evaluate the combination 

^.V^V) = (foo + f,,) = Too + (2) 

where we have choosen the null vector in z direction: = (1, 0, . . . , 0, 1). In 
a next step we expand the scalar field $ in terms of momentum modes tpn- These 
momentum modes are defined as eigenmodes of the Laplacian in the presence of 

(-V^ + a(x)) iJn{x) - kltPn(.x), (3) 

with fc^ denoting the eigenvalue. We now introduce the propagator G{x, x' , k) of 
the corresponding Helmholtz equation 

(-V^ + (j{x) - k^) G{x, x', k) = 5{x - x') (4) 

which has the spectral representation 

G(.y.k)^Y.f^,^ (5) 

n ^ 

The components of the energy-momentum tensor in Eq. ([2| can now be expressed 
solely in terms of G{x,x' ,k). As we are interested in the effect of the background 
potential (t{x) in comparison to the vacuum case cr = 0, we normalize the EMT with 
respect to the trivial vacuum by subtracting the free Green's function Gq[x,x' ,k) 
obtained from Eq. Q for a = 0. This corresponds to removing the cr-independent 
divergence, 

f°° dk / I \ 
{Tooix,t))= lim / —ik^ + -V-iV + V)]lm[iG-Go)ix,x',k)], (6) 







TT V 2 



dk / 1 

{T,,{x,t))^ lim / — a,a,.--V-(V + V'))lm[(G-Go)(a;,a;',A:)]. (7) 







2 



The whole formalism works, of course, for any field configuration 17(0;). In the fol- 
lowing, we focus on a backgrounds that induce Dirichlet BCs on dV. For this, we 
schematically use 

cr{x) ^ XS {x - XQv) , XQ-D&dV, (8) 

and take the Dirichlet hmit A ^ 00 ^o. The EMT in Eq. Q and Q is then finite 
on V where the potential a vanishes. There remain, however, divergences on the 
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boundary dT> due to the Dirichlet limit. These can be interpreted as the infinite 
amount of energy needed to constrain a scalar quantum field on all momentum 
scales to satisfy the Dirichlet constraint on the boundarjCHl, 

Interpreting the eigenvalue problem (|3]) as a ficticious quantum mechanical 
Schrodinger problem with Hamiltonian H = — + cr(a;), the Green's function 
corresponds to a quantum mechanical propagator Fourier transformed to energy 
space, which can be written as a Feynman path integral in position space. The 
vacuum-subtracted path integral representation of the propagator then results in 

oo x=x{T) = .2 / \ 

iG-Go){x,x',k)=i dse'^'' / Vx{T)eo Me o _i 



x'=x{0) 

dTe-^'^- / Vxir)e » Me ' -1 , (9) 



x' 

where s is a ficticious quantum mechanical time, so-called propertinie. We have 
formally made Wick rotations both in the s and k planes so that s = —iT and 
— ik to cast the expression above in Euclidean form. 

In contrast to previous calculations of effective interaction energies for the 
Casimir effect and similar setups which required path integrals only over closed 
loops, here we need to compute path integrals of open worldlines from x' to x. 
This can conveniently be done with the d loop algorithnpH. The path integral is 
normalized here implicitly such that the free path integral yields the standard free 
propagator, 

^ T , i2 (x-x')^ 

Vx{T)e = -. (10) 

This normalization is made explicit by the following shorthand notation 

e 4T I -j dr^ 

(O),,, := / Vx{r)e « O. (11) 

Then Eq. ^ can be written as 



(x-x')2 / T 

1.2 e ' IT / -JdTa{x) 



{G-Go){x,x\k = -tkE)= I dTe'^''^- — { e '> ^'-l) ■ (12) 

(AttT) '■ 







-/drff(x) \ _ J -1 if 2;(r) intersects a boundary 1 _ 



In the Dirichlet limit the last term can be simplified 

— 1 if x(t') intersects a bour 
otherwise J 

(13) 

which means that only loops x{t) that violate the boundary conditions lead to 
deviations from the trivial vacuum and thus contribute to the path integral. The 
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quantity Tmin denotes the minimum propertime needed for a given worldline to 
propagate from x' to x and simultaneously intersect a boundary in between. If T 
is too small the particle has to propagate on a straight line from a;' to a;. If T 
is sufficiently long the diffusive Brownian motion process behind the path integral 
creates random detours that can eventually intersect the boundaries. We may choose 
our integration contour in the kE plane in such a manner that the imaginary part 
of the integral is already taken into account by carrying out the integrations 
in ([6| and Q along the real axis, indeed this is the effect of the Wick rotation. 
These considerations lead us to a compact representation of the EMT components 
required, e.g., for the null energy condition: 

CO 

(foo(a;,t))(/) ^ — 4^ / -^e-T (9 (T - r„,i„)),,, , (14) 
(47r)^ J 2T^ 

CO 

^ 1 /" dT^ ix — x'^"^ 

{TUx:t)){I) = - —^dA- / -^e-^^ {Q{T~T^^)),^, , (15) 

(47r)T" J 
^ ' 

oo 

(foo(x,0)(//)= ---4TrV.(V + V') / ^e-^(e(T-T^i„)),,, 
(47r)~ J 2T— 

= -(f,,(a;,i))(//), (16) 

where the limit a; — > is implicitly understood from now on. From the last identity, 
it is, for instance, obvious that the NEC only involves ^T^,yV^V'''^ = TooiI)+T^ziI). 
We emphasize that the angle brackets are used in two different contexts: on the l.h.s. 
(...) denotes the vacuum expectation of a composite quantum operator, whereas 
on the r.h.s. (. . . )^^, denotes the expectation value of the path integral, that is the 
average over an ensemble of open worldlincs. 



3. Energy momentum tensor for a single plate 

With these general considerations, we can now compute the EMT components in 
the case where the boundary is one single plate at z = with the z axis being its 
normal. 



3.1. Analytical calculation for the single plate 

Analytical results for the single plate are easily obtained through explicit calculation 
of the Green's function using the method of images. After regularization one finds 

{G~Go){x,x';k)^~ip{^,x') (17) 

in which ip{S^,x') denotes the field at x' :— (x'|,z') produced by a point source 
sitting at :— (xy,— z). Also x := (a;||,z), with z and z' being the distances from 
the plate. The explicit forms of the field in two and three dimensions are given 
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by (pd=2 = ^Koi-ik\x - x'\) and (pd^s = lt\^-V\ , respectively. Using them it is 



straightforward to show that 



^oo)^^^(I)= ~, (18) 

1 1 

d=2''' 16^^' \''-Vd=3''' 327r2z4 



T..). JI) =-T7^Z^, {T..), JI) (20) 



1 


1 


327r2 




3 


1 


327r2 




3 


1 



3.2. WorZdZme calculation for the single plate 



For the computation of the expectation value in Eq. (14) and (16 1, we need to 
evaluate the worldline Q function, i.e., we need to find the condition under which a 
loop in our ensemble intersects the plate (see Fig. [l] for a sketch) . 




Fig. 1. Sketch of the single plate setup. 



After rescaling the loops for numerical convenience' 

':{Tt) ^ Vfy (t), t e [0,1] / Vx{T)e ° ^ / Vy(t)e ° 

Jx' Jx'/VT 

(21) 

we find in the limit x' — >■ x, that is, z' — z, the intersection condition 

Vfy^ + z<o =^ e (T - r„,i„) = e(T-^y (22) 

where z denotes the distance from the plate and y_ denotes the point of a loop 
that is closest to the plate, see Fig. [T] The quantity y_ is the only information we 
need about each loop. Since the worldline distributions factorize with respect to 
their position space components, we only need 1-dimensional loops. In d = 2 the 
Too component of the EMT then reads 

(Too) ^ (Too) (I) + (Too) (II) ^ - i\y-\)] (23) 
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Table 1. Worldline results for the single plate in d = 2. 



ppl, N 


(|y-l> 


2 


(|y-l'> 


3^ 
4 


2", 10* 


0.8789 ± 0.0046 


0.8862 


1.3018 ±0.0210 


1.3293 


2i*, 5 ■ 10^ 


0.8841 ± 0.0007 


0.8862 


1.3236 ±0.0029 


1.3293 



with the corresponding analytical values in the second line. In Table [T] we compare 
our analytical and numerical results for two different ensembles with N loops per 
ensemble and ppl points per loop. The error corresponds to the statistical uncer- 
tainty only and it decreases like 1/VN as expected. The systematic error that is 
not displayed explicitly is governed by the discretization of the loops (i.e. by ppl). 
A larger number of points per loop improves our results. As is known from previous 
calculations of the Casimir energy using worldline numerics the numerical values 
approach the analytical values from below, that is the discretized loops are smaller 
than the continuous ones. 

The same analysis can be done for three spatial dimensions where we find 




Too) = j^i,\^^-^-^^] (25) 

(26) 

The structure of the result is, of course, identical but the exponents of our loop 
variable and coefficients have changed. In Table[2]we show our worldline results. Just 
as in the two-dimensional case we again see the decrease of errors with increasing 
number of loops and ppl. In both cases, these results confirm the violation of the 
weak energy condition by the single-plate configuration in agreement with Ref. |31 



Table 2. 


Worldline results for the 


single plate in d = 


3. 


ppl, N 




1 




2 


2i4, 10* 


0.9841 ± 0.0100 


1.0000 


1.9610 ±0.0488 


2.0000 


218, 5 ■ 10^ 


0.9964 ± 0.0014 


1.0000 


1.9901 ± 0.0063 


2.0000 



As a further check of our algorithm we compute the null energy condition along 
the z axis 

which is also violatecP. Comparing with the analytical results again shows that our 
algorithm reproduces the expected values within the numerical precision set by the 
parameters of our loop ensemble, viz. the number of points per loop ppl and the 
number of loops N. 
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4. Energy density for two parallel plates 



Let us now turn to the classic Casimir configuration of two infinite parallel plates 
separated by a distance a imposing Dirichlet BCs on the fluctuating scalar field. 
However, because our energy-momentum tensor is finite everywhere except for dT> 
we are not only able to compute energy conditions or EMT components but actually 



also parts of components. In the following, we focus on Too(/), cf. Eq. (14) 



4.1. Analytical calculation for parallel plates 

As before Too(/) may be calculated by means of the traditional Green's function 
approach and the method of images which gives 

(G-Go)(x,x';fc)= ^ ^{x„x')-Y,^i^,^^')- (28) 

9ez/{o} ?ez 

Here the vectors Xg i^w^z + 2qa), :— (a;||,— z + {2q + l)a) were introduced. 
This time, z measures the distance from the plates in a coordinate system centered 
in the middle of the plates. Introducing the dimensionless coordinate C = z/a and 



inserting the expressions for ip given above Eq. ( 18 ) yields for the region between 
the plates where |C| < 1/2 

with Cr,h denoting the Riemann and Hurwitz zeta functions, respectively. 



4.2. Worldline calculation for parallel plates 

The numerical worldline computation of 

CO 

Too) (I) = 7-4+1 / r5±F (0 (T - Ilnin)) (30) 

^ ' 

for the Casimir plates is similar to the single-plate case. However, as there are 
now two plates we need two intersection conditions as well. Accordingly we must 
compute two points of each loop, viz. the points y± which are closest to the plates, 
see Fig. [2] The corresponding value of Tmin is then easily found to be 



a mm 



(31) 



The analytical findings of Eq. ( 29 ) compare favorably with the numerical data 
displayed in Fig. [3j 
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Fig. 2. Sketch of Casimir's parallel-plates configuration in units of a = 1. 



{T-oo)(l) 




C 

Fig. 3. EMT component Too(/) in between two Casimir plates at a distance a in two and three 
spatial dimensions, C = zjeL. The analytical result Eq. ( |29| (solid lines) is very well approximated 
by the numerical data. 

4.3. Conclusion 

We have studied Casimir energy-momentum tensors induced by a fluctuating mini- 
mally coupled scalar field obeying Dirichlet boundary conditions. We have demon- 
strated that this general problem can be formulated with the aid of the worldline 
approach to the Casimir effect. Contrary to simple interaction energy computations, 
the worldline approach for such composite operators now has to deal with open 
worldlines arising from the point-splitting procedure. As a proof of principle, we 
have been able to show that the numerical worldline approach to composite opera- 
tors is able to reproduce analytically known results for the single- and parallel-plate 
cases at low cost and high efficiency. In particular, the single-plate case already 
shows violations of the weak as well as the null energy condition. Of course, the 
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main goal of this project is to study more general classes of Casimir configurations 
which particularly allow for non-trivial tests of the averaged null energy condition. 
Work in this direction is in progress. 
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